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(54) Channel estimation using a sliding window technique 

(57) The invention provides a method for channel 
estimation in mobile radio communications which adap- 
tively compensates for channel distortion on a block-by- 
block basis. The discrete-time channel impulse 
response is initially estimated with a given length and 
then truncated by using a sliding window. A cost func- 
tion associated with the window is measured as the 
length and position of the window is adjusted over the 
channel impulse response and the cost function is com- 
pared with a threshold. The invention provides means to 
and a method for adaptively adjusting the length of the 
window L T and the corresponding number of states in 
the equalizer, 2 LT - 1 if appropriate. 
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Description 



[0001] This invention relates to channel estimation methods and apparatus in mobile radio communications and in 
particular to a receiver which adaptively compensates for channel distortion on a block-by-block basis. 
5 [0002] In digital mobile radio communications, transmission channels suffer from severe distortion due to frequency 
selective fading. In addition, channel characteristics are normally time-varying due to the relative motion of fixed and 
mobile stations. Therefore, in order to allow for reliable transmission, the receiver must be able to estimate and com- 
pensate for channel distortion on a block-by-block basis. 

[0003] Various channel estimation and channel equalization methods have been proposed in literature and are com- 
10 monly used in practical systems such as mobile cellular communication systems employing the European wireless dig- 
ital cellular standard "GSM". In most cases the receiver performs channel equalization on the received signal using 
Maximum Likelihood (ML) or Maximum A Posteriori (MAP) probability data estimation, based on the knowledge of the 
Channel Impulse Response (CIR). Most practical systems employ training sequences to enable the CIR to be esti- 
mated before the equalizer star-up. Fast time varying, fading channels require the changing channel response to be 
is tracked and adjusted dynamically by the receiver for the duration of the received signal. Tracking of the CIR may be per- 
formed by means of decision directed algorithms, where tentative decisions from the equalizer are used to update the 
initial CIR estimate. Examples of receiver systems which perform channel estimation and channel equalization may be 
found in the following articles: "Bit Synchronisation and Timing Sensitivity in Adaptive Viterbi Equalizers for Narrowband 
TDMA Digital Mobile Radio Systems", A. Baier, G. Heinrich and U. Wellens. Proc. IEEE Vehicular Technology Confer- 
so ence. June 1988, pp 377-384 [Reference 1]; "Correlative and Iterative Channel Estimation in Adaptive Viterbi Equalize 
ers for TDMA Mobile Radio", JTG-Fachbericht No. 107, VDE Verlag, April 1989, pp 363-368 [Reference 2]; "Simulation 
and Hardware Implementation of a Viterbi Equalizer for the GSM TDMA Digital Mobile Radio System", A. Baier, G. 
Heinrich, P. Shoeffel and W. Stahl, Proc. 3 rd Nordic Seminar on Digital Land Mobile Radio Communications, September 
1988 pp 13.7.1 - 13.7.5, [Reference 3]. 
25 [0004] An equalizer of a given complexity can only cope with a certain delay spread of signal paths. Therefore, the 
estimated CIR is usually truncated and this truncated CIR estimate is used in a reduced state equalizer, of the type 
described in "Mobile Radio Communications" by R. Steels, pub. By Rentech Press and IEEE Press, pp 560-575. The 
reduced state equalizer selects an appropriate segment of the estimated CIR by sliding a window of length L T over the 
whole estimated response, calculating the energy contained within the window at each window position and identifying 
30 the window position where the energy within the window is maximum. 

[0005] However, if the actual CIR is either longer or shorter than L T samples, the above approach results in a degra- 
dation of the equalizer performance. In fact as it is intuitive, assuming a channel response shorter than required pro- 
duces a residual ISI which will affect the data estimation. On the other hand, assuming an impulse response longer than 
necessary means making use of one or more taps of the estimates channel which are essentially noise. There is thus 
35 a requirement for a reduced-state equalizer which provides improved equalizer performance. 

[0006] According to a first aspect of the invention there is provided a method of estimating channel impulse response 
in a communication system comprising: 

initially estimating the discrete-time channel impulse response with a given discrete length; and 
40 truncating the channel impulse response by using a sliding window, and characterised in that a cost function asso- 
ciated with the window is measured as the length and position of the window is adjusted over the channel impulse 
response and the said cost function is compared with a threshold. 

[0007] The overall performance of the equalizer is improved by neglecting the portion of the estimated CIR which is 
45 likely to correspond to just estimation noise. 

[0008] The invention provides means to and a method for adaptively adjusting the length of the window L T and the 
corresponding number of states in the equalizer, 



if appropriate. This may be achieved by computing a cost function for different lengths and positions of the window, for 
example by measuring the out-of-window power for each window, and comparing that cost with a suitable threshold. 
[0009] Once L samples CIR have been obtained by standard means, the position of the truncated CIR of different 
55 lengths L T is then obtained by performing several maximum energy searches based on sliding windows of length L T . 
For each window length, the out-of-window energy of the estimated CIR is compared with a suitable threshold. The 
threshold may be fixed, or it may be obtained from a previous estimate of the SN R (based on the non-truncated channel 
estimate). The invention may also be implemented by comparing the energy of the individual CIR taps with a suitable 
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threshold. Again, the threshold may be fixed, or it may be obtained from a previous estimate of the SNR (based on the 
non-truncated channel estimate). 

[001 0] A detailed description of a practical digital radio receiver is described below, by way of example, and with ref- 
erence to the following figures in which: 

5 

Figure 1 shows in outline a typical GSM digital radio receiver; 
Figure 2 illustrates the GSM "normal" burst format; 
Figure 3 illustrates a sliding window; 

10 [001 1 ] A typical implementation of a digital radio receiver is shown in Figure 1 and comprises an equalizer 1 0, dein- 
ter leaver 15 and a channel decoder 16. The equalizer includes a channel estimator 11 and data estimator 14. The 
channel estimator includes an initial channel estimator 12 and a channel truncator 13. 
[001 2] The discrete-time received signal as received by the equalizer can be written as 

15 

r(kh^b(k-l)h(l) + n(k) (1) 



20 where b(k) e {-1 ,1} are the transmitted data symbols or the (known) training sequence symbols, h(l) I = 1 , 2 L rep- 
resents the taps of the Channel Impulse Response (CIR) and n(k) indicates white Gaussian noise with zero mean and 
variance a 2 . 

[0013] The equalizer must first estimate (12) the CIR h(l), before beginning the data estimation process (14). In some 
cases, e.g. in a GSM standard receiver, the initial CIR estimation is commonly performed by means of correlative chan- 
25 nel sounding, see for example the above mentioned references 1 and 3. The taps of the CIR estimate are obtained by 
correlating the received signal r(k) with N= 16 bits b(k) out of the 26 bits of training sequence, shown in figure 2, and 
represented by: 

h(l)=(VN)xp{i)r{l + i) (2) 

30 

where 

/»0. f,.„ N-1; 

35 [0014] The initial CIR estimate can generally be performed by ML channel estimation, giving: 

A = [A(0),A(l),...,A(/:-l)] r =(B'B)- l B r r (3) 

40 

where 

r = [r(0),r(l),...,r(AT-l)] r 
B = [b(0),b(l),...b(tf-l)] r 

50 

b(Q-[£>(Q.f>(M) b(i-L+l)) T 

[001 5] It can be seen that, due to the good autocorrelation properties of the GSM training sequence 

55 

8*8 = NI 
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and equation 2 is the particular case of the more genera! ML technique (equation 3). Once L samples (h(l)) of the CIR 
are estimated, an appropriate segment of samples (h for L T < L) is selected (in the channel truncator 1 3) by a max- 
imum energy search based on a sliding window at length L T , see reference 1, and as shown in figure 3. This approach 
produces a bit synchronisation with a resolution equal to the sampling period T and results in a reduced complexity 
5 equalizer with a trellis of 
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[001 6] Assuming that for the training sequence bits b(k) 



is substituting (1) in (2) yields 
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h(l)=h(l)+v(l). 1=0,1 L-1. (5> 



where 



AM 

*M" rZ WW) '=0. -AM , (6) 



[0017] Suppose that Lj taps of h(1) are selected as a CIR estimate hj{l) t when the actual channel response is 
shorter than L T bit intervals. From (5), it is clear that if h(l) = 0 for some index I, the corresponding sample h is just 
noise, and its use can only degrade the equalizer performance. An improvement may be achieved by adaptively adjust- 
ing the parameter L T on a block-by-block basis, and consequently adjusting the number of states 



2 v 



of the equalizer, in order to match the actual length of the channel impulse response. 
35 [0018] A suitable cost function, such as the mean square error of the truncated channel estimate, may be used to 
estimate length of the actual CIR and is given by 

J('o.^>=£{lW0-M0l 2 }. (7) 

40 /=0 

where h T (l) = h(l) for / e W = [/ 0 ,/ 0 + L T - 1] , and h T (l) = 0 for / g W. Taking into account (5), he above quantity 
can be written as 



J(/ 0 .L T )o £ E{|v(0l 2 }+ 2 E[\h(l)\ 2 } (8) 

Therefore, choosing the value of (/ 0 , L T ) which minimises the cost (7) implies a trade-off between reducing the estima- 
tion noise and representing the relevant IS! terms. From (5) and (7) one finally obtains 
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J(/ 0 ,Z. T ) = £ E{|v(/)| 2 }+ 2 |*</)| 2 - 2 E{k(/)| 2 } (9) 

teW ** w l*W 

= (2L r -L)E{| V (/)| 2 }+ S IW)! 2 

[001 9] One can observe that, for a fixed window length L T , minimising (9) is equivalent to the maximum energy search 
described above. 

10 [0020] However, in comparing the cost associated with different window lengths Lj, the above approach requires an 
accurate estimate of the channel noise variance a 2 , to be employed to evaluate the mean squared value of the estima- 
tion error v(k) as defined in (6). From this, as an alternative to (9), one can consider the following process. Assume that, 
for each window length L Tf the best window position / 0 is determined by a maximum energy search. Once a value of /q 
is associated with each length L T , the best L T can be chosen by comparing the out-of-window energy 
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5>{|v(/)i 2 



20 corresponding to the different values of L T with a suitable threshold. In the case where the out-of-window taps of the 
estimated channel do not contain any relevant ISI terms, using (5) and (6) and by denoting with E L and E LT the energy 
of h(f) and hrfl) respectively, for signal-to-noise ratios 

25 ^A/ 0 = (Zl /? ( / )| 2 )/a 2 >Y 



one has 

30 

E L -E Lr = (L-L 7 )(a 2 //V) < (L-L T ) (E LT /yN). (10) 



[0021 ] Therefore, for each burst, the quantity E L - E LT is computed and the validity of (6) for different values of L T is 
35 checked. The lowest L T satisfying the inequality (6) is an estimate of the length of the actual channel response, and may 
be then used both as a length of the CIR and for setting the number of states in the equalizer. 
[0022] The effectiveness of the invention has been asses by computer simulation for the case of a GSM receiver. The 
GMSK transmitted data symbols are obtained from the source bits by a rate 1 :2 convolutionally encoding and interleav- 
ing , according to the GSM specifications for the TCH/FS class lb bits. At the receiver, the channel estimation is per- 
40 formed by using N=16 bits out of the 26-bits training sequence midamble of the GSM normal burst, according to (2). 
The original L = 8 taps estimate are then truncated to L T £ 6 samples according to the rule (9) or (10). For data estima- 
tion, we employ a symbol-by-symbol Max-Log-MAP equalizer on the ISI trellis with state complexity 



45 



The equalizer soft-output data are deinterleaved and then fed to a symbol-by-symbol-Max-Log-MAP convolution chan- 
nel decoder. 

[0023] Two channel profiles from "Digital cellular telecommunications system (phase 2+): Radio transmission and 
so reception." ETSI GSM 05.05 (version 5.2.0), July 1996 are considered. The first is the GSM typical urban area profile 
(TU50) with a relative speed of 50 Km/h and a delay spread of about S^sec. The second propagation condition is the 
GSM typically hilly terrain profile (HT100), with a relative speed of 100Km/h and a delay spread of about 20jisec. Con- 
sidering the GSM symbol interval T = 3.69 jisec, and taking into account that the GSMK pulse duration is about 3T, the 
above propagation conditions correspond to the two opposite cases where the invention provides a significant (TU con- 
55 ditions) and a negligible (HT conditions) performance improvement respectively. 

[0024] To measure the performance of the different channel estimation algorithms, the means square error between 
the true and estimated channel responses are considered. In the case of the TU50 conditions, both strategies (9) and 
(10 produce a gain at the equalizer output and at the channel decoder output. For the HT100 profile, no significant 
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improvement is achieved and the BER is similar to a conventional receiver. 
Claims 

5 1 . A method of estimating channel impulse response in a communication system comprising: 

initially estimating the discrete-time channel impulse response with a given length; and 
truncating the channel impulse response by using a sliding window, and characterised in that a cost function 
associated with the window is measured as the length and position of the window is adjusted over the channel 
10 impulse response and the said cost function is compared with a threshold. 

2. A method as claimed in claim 1 characterised in that the cost function is computed for each different length and/or 
position of the window. 

is 3. A method as claimed in claim 1 characterised in that the out-of-window power is used as the cost function. 

4. A method as claimed in claim 1 characterised in that the length and position of the window given the best result, 
when compared with the threshold, is chosen as the channel impulse response estimate. 
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